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o ,■ Abstract 

In the present work we study both thermal and non-thermal leptogenesis in 
all neutrino mass models describing the presently available neutrino mass 
patterns. We consider the Majorana CP violating phases coming from right- 
handed Majorana mass matrices to estimate the baryon asymmetry of the 
universe, for different neutrino mass models namely degenerate, inverted hi- 
CS1 ■ erarchical and normal hierarchical models, with tribimaximal mixings. Con- 

^ . sidering two possible diagonal forms of Dirac neutrino mass matrix as either 

charged lepton or up-quark mass matrix, the right-handed Majorana mass 
matrices are constructed from the light neutrino mass matrix through the 
inverse seesaw formula. Only the normal hierarchical model leads to the best 
predictions for baryon asymmetry of the universe, consistent with observa- 
tions in both thermal and non-thermal leptogenesis scenario. The analysis 
though phenomenological may serve as an additional information in the dis- 
crimination among the presently available neutrino mass models. 
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1 Introduction 

The existence of heavy right-handed Majorana neutrinos in some of the 
left-right symmetric GUT models, not only gives small but non vanishing 
neutrino masses through the celebrated seesaw mechanism [1], it also plays 
an important role in explaining the baryon asymmetry of the universe [2] 
Yb ={6-1-o%) 1CT 10 . Such an asymmetry can be dynamically generated if 
the particle interaction rate and the expansion rate of the universe satisfy 
Sakharov's three famous conditions [3]. Majorana right-handed neutrinos 
satisfy the second condition i.e, C and CP violation as they can have an 
asymmetric decay to leptons and Higgs particles, and the process occurs at 
different rates for particles and antiparticles. The lepton asymmetry is then 
partially converted to baryon asymmetry through the non-perturbative elec- 
troweak sphaleron effects [4,5]. In such thermal leptogenesis the right-handed 
neutrinos can be generated thermally after inflation, if their masses are com- 
parable to or below the reheating temperature Mi < Tr. This allows high 
scale reheating temperature Tr > 10 9 GeV[ 6]. In non-thermal leptogene- 
sis [7] it is possible to produce lepton asymmetry by using the low reheating 
temperature, where the right-handed neutrinos are produced through the di- 
rect non-thermal decays of the inflaton. This is particularly important for 
supersymmetric models where gravitino problem [8] can be avoided provided 
the reheating temperature after inflation is bounded from above in a certain 
way, namely T R < (10 6 - 10 7 ) GeV. 

In order to calculate the baryon asymmetry from a given neutrino mass 
model, one usually starts with the light neutrino mass matrices ttill and then 
relates it with the heavy Majorana neutrinos M^ and the Dirac neutrino 
mass matrix tulr through inverse seesaw mechanism in an elegant way. We 
consider the Dirac neutrino mass matrix tjilr as either the charged lepton 
mass matrix or up quark mass matrix for phenomenological analysis. The 
complex CP violating phases are usually derived from the MNS leptonic mix- 
ing matrix. In the present work we are interested to consider the complex 
Majorana phases which are derived from the right-handed Majorana mass 
matrix Mrr, in the estimation of baryon asymmetry of the universe. We 
wish to consider the left-handed light Majorana neutrino mass matrices tull 
which obey the \x — r symmetry[9] where tribimaximal mixings[10] are re- 
alised, for all possible patterns of neutrino masses, viz, degenerate, inverted 
hierarchical and normal hierarchical mass patterns. We first parametrise the 
light left-handed Majorana neutrino mass matrices which are subjected to 



correct predictions of neutrino mass parameters and mixing angles. The cal- 
culation of baryon asymmetry of the universe in the light of thermal as well 
as non-thermal leptogenesis, may serve as an additional information to fur- 
ther discriminate the correct pattern of neutrino mass models and also shed 
light on the structure of Dirac neutrino mass matrix. 

In section 2, we briefly mention the formalism for estimating the lep- 
ton asymmetry in thermal leptogenesis through out of the equilibrium decay 
of the heavy right-handed Majorana neutrinos, followed by numerical cal- 
culation and results. Section 3 is devoted to non-thermal leptogenesis and 
numerical predictions. Finally in section 4 we conclude with a summary and 
discussion. Important expressions related to tull which obey \x — r symmetry 
for three neutrino mass models, are relegated to Appendix A. 

2 Baryon asymmetry of the universe in ther- 
mal leptogenesis 

The canonical seesaw formula (known as type-I)[l] relates the left-handed 
Majorana neutrino mass matrix ijill and heavy right handed Majorana mass 
matrix Mrr in a simple way 

m LL = -m LR Mg l R m T LR (1) 

where tulr is the Dirac neutrino mass matrix. For our calculation of lepton 
asymmetry, we consider the model[5,ll] where the asymmetric decay of the 
lightest of the heavy right-handed Majorana neutrinos, is assumed. The 
physical Majorana neutrino Nr decays into two modes: 

N R - l L + 0t 

where II is the lepton and II is the antilepton and the branching ratio for 
these two decay modes is likely to be different. The CP-asymmetry which 
is caused by the intereference of tree level with one-loop corrections for the 
decays of lightest of heavy right-handed Majorana neutrino Ni, is defined by 

[5,12] 

e = ^£ 
r+r 



where T = T(Ni — > Il4>) and T = T(Ni — > /^0) are the decay rates. A 
perturbative calculation from the interference between tree level and vertex 
plus self energy diagrams, gives [13] the lepton asymmetry e\ for non-SUSY 

CclSG clS 
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where /(x) and ^(x) represent the contributions from vertex and self-energy 
corrections respectively, 

f(x) = Vx[-1 + (x + l)/n(l + -)], 

x 



'X 

g(x) 



x — 1 

For hierarchical right-handed neutrino masses where x is large, we have the 
approximation [2], f(x) + g{x) ~ 27s- This simplifies to 



ei 



167T 
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where /i = ttilr/v is the Yukawa coupling of the Dirac neutrino mass matrix 
in the diagonal basis of Mrr. In term of light Majorana neutrino mass matrix 
will, the above expression can be simplified to 

For quasi-degenerate spectrum i.e., for Mi ~ M 2 < M 3 the asymmetry is 
largely enhanced by a resonance factor and in such situation, the lepton 
asymmetry is modified [14] to 

„ 1 MM , (4) 



where 



"" ~ (Mpi^P+rfMf allCL X 2 ~ 8^ 



It can be noted that in case of SUSY, the functions f(x) and g(x) are given 
by f(x) = yfxln(l + -) and g(x) = -^|; and for large x one can have 
f(x) + g(x) ~ -t=. Therefore the factor | will appear in place of y| in the 
expression of CP asymmetry [2]. 

The CP asymmetry parameter e\ is related to the leptonic asymmetry 
parameter through Yl as 

_ 3 

s i 9*i 

where ul is the lepton number density, ni in the anti-lepton number density, 
s is the entropy density, «$ is the dilution factor for the CP asymmetry e^, 
and g*i is the effective number of degrees of freedom at temperature T = Mj. 
The baryon asymmetry tib produced through the sphaleron transmutation 
of Yi, while the quantum number B — L remains conserved, is given by [15] 

n ~f - cr-t - <£t* (6) 

where 

= 8N F + AN H 

22N F + 13N H ' { ' 

Here iVp is the number of fermion families and iV# is the number of Higgs 
doublets. Since s = 7.04n 7 the baryon number density over photon num- 
ber density n 7 corresponds to the observed baryon asymmetry of the Uni- 
verse [16], 

Y sm = { V±)SM „ dKi€i (g) 

n 7 

where d ~ 0.98 x 1CT 2 is used in the present calculation. In case of MSSM, 
there is no major numerical change with respect to the non-supersymmetric 
case in the estimation of baryon asymmetry. One expects approximate en- 
hancement factor of about v / 2(2y / 2) for strong (weak) washout regime[2]. 

In the expression for baryon-to-photon ratio K\ descrides the washout 
of the lepton asymmetry due to various lepton number violating processes. 
This efficiency factor (also known as dilution factor) mainly depends on the 
effective neutrino mass rh\ 

" Ll ~ Mi 



where v is the electroweak vev, v = 174GeV. For 10 2 eV < rh\ < lO 3 ^, 
the washout factor K\ can be well approximated by[12,17] 



^i(^i) = 0.3 
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Vfl\ 



log 



mi ' a6 
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(9) 



We adopt a single expression for K\ valid only for the given range of mi[17,18,19]. 

2.1 Numerical calculations and results 

To compute the numerical results, we first choose the light left-handed Majo- 
rana neurino mass matrix tjill proposed in Appendix A. These mass matrices 
obey the \x — r symmetry[9] which guarantees the tribimaximal mixings[10]. 
The input parameters are fixed at the stage of predictions of neutrino mass 
parameters and mixings given in Table 1. These results are consistent with 
the recent data on neutrino oscillations. 

For the calculation of baryon asymmetry, we then translate these mass 
matrices to M RR via inversion of the seesaw formula, M RR = —rn^ R rril\rn LR . 
We choose a basis U R where M RR = U R M RR U R =diag(Mi, M 2 , M 3 ) with real 
and positive eigenvalues. We then transform diagonal form of Dirac mass 
matix, mL R =diag(\ m , \ n ,l)v to the U R basis: ttilr — > m' LR = rriLRU R Q 
where Q = diag(l,e ta ,e tf3 ) is the complex matrix containing CP-violating 
Majorana phases derived from M RR . Here A is the Wolfeinstein paramater 
and the choice (m, n) in m^ R gives the type of Dirac mass matrix. For ex- 
ample, (6,2) for charged-lepton type mass matrix and (8,4) for up-quark 
type mass matrix. In this prime basis the Dirac neutrino Yukawa cou- 
pling becomes h! = -^^ which enters in the expression of CP-asymmetry 
e±. The Yukawa coupling matrix hi also becomes complex, and hence the 
term Im(h'h)ij appearing in lepton asymmetry e\ gives a non-zero contribu- 
tion. A straightforward simplification shows that (h'h)fj = (QI^Q^^ + 
{Qii) 2 Q33^-3 where i?2,3 are real parameters. After inserting the values 
of phases the above expression leads to Im{Wh)\j = — [i?2sin2(a — /3) + 
i?3 sin 2a] which imparts non-zero CP asymmetry for particular choice of 
(a,/3). 

In our numerical estimation of lepton asymmetry, we choose some arbi- 
trary values of a and (3 other than tt/2 and 0. For example, light neutrino 
masses (mi, — m 2 , m 3 ) lead to M RR = diag(M 1 , — M 2 , M 3 ), and we thus fix 
the Majorana phase Q = diag(l,e (ia \e ii/3) ) = diag(l,e i(7T/2+7v/4 \e i7r/4 ) for 



Type 


Am 2 21 [10- ! 


W 2 ] 


Am| 3 [10- ; 


'eV' 2 } 


tan 2 # 12 


sin 2 2# 23 


sin 6>i3 


Deg.(IA) 
Deg.(IB) 
Deg.(IC) 


7.8 
7.9 
7.9 




2.6 

2.5 
2.5 




0.5 
0.5 
0.5 


1.0 
1.0 
1.0 


0.0 
0.0 
0.0 


Inh.(IIA) 
Inh.(IIB) 


7.3 

8.5 




2.5 
2.3 




0.5 
0.5 


1.0 
1.0 


0.0 
0.0 


Nh.(IIIA) 
Nh.(IIIB) 


7.1 
7.5 




2.1 
2.4 




0.5 
0.5 


1.0 
1.0 


0.0 
0.0 



Table 1: Predicted values of the solar and atmospheric neutrino mass-squared differences 
for tan 2 #12=0.50, using m^i given in the Appendix A. 

a = {tt/4: + 7r/2) and (3 = ir/4. The extra phase 7r/2 in a absorbs the 
negative sign before heavy Majorana mass M 2 . In our search programme 
such choice of the phases leads to highest numerical estimations of lepton 
CP asymmetry. 

In Table 1 we give the predictions on Am 2 ,! and Am 2 ^ of these seven 
neutrino mass models under consideration in Appendix A. They obey ft — r 
symmetry and predict tribimaximal mixings in addition. In Table 2 the three 
heavy right-handed neutrino masses are extracted from the right-handed Ma- 
jorana mass matrices so constructed through inverse seesaw formula, for three 
choices of diagonal Dirac neutrino mass matrices. We get degenerate spec- 
trum of heavy Majorana masses for normal hierarchical model and this allows 
us to use resonant leptogenesis formula. The corresponding baryon asymme- 
try Yb are estimated in Table 3 and this shows that only normal hierarchical 
model predict reasonable values whereas inverted hierarchical model (IIB) 
nearly misses the observational bound. Degenerate models predict too low 
baryon asymmetry. 

Our estimatated baryon asymmetry for normal hierarchical model (IIIA. 
IIIB) lies between 9.27 x 10~ 9 with Dirac neutrino mass matrix as charged 
lepton mass matrix (6,2), and 7.28 x 10 -11 for the up-quark mass matrix 
(8,4). This hints a possible choice of Dirac neutrino mass matrix lying be- 
tween these two e.g., mm = diag.(X s ,X 2 ,l)v. As emphasised earlier, our 
starting point is the neutrino mass matrix which satisfies the observed neu- 
trino mass parameters and mixings. The values of input parameters are fixed 
at this level before applying to the calculation of baryon asymmetry. The 
whole calculation is performed in a consistent way. 



Type 


m,r 


Mi 


M 2 


M 3 


IA 


'6,2 N 


1.22 x 10 8 


-6.01 x 10 11 


2.59 x 10 13 


IA 


'8,4^ 


9.86 x 10 5 


-5.03 x 10 9 


2.51 x 10 13 


IB 


'6,2 N 


4.05 x 10 7 


6.16 x 10 11 


7.60 x 10 13 


IB 


'8,4^ 


3.28 x 10 5 


4.99 x 10 9 


7.60 x 10 13 


IC 


'6,2^ 


4.05 x 10 7 


-6.69 x 10 12 


6.99 x 10 12 


IC 


'm: 


3.28 x 10 5 


-4.83 x 10 11 


7.84 x 10 11 


IIA 


'6? 


3.29x 10 8 


9.73X10 12 


6.25xl0 16 


IIA 


'8,4^ 


2.63x 10 6 


7.94xl0 10 


6.21xl0 16 


IIB 


'6,2^ 


-9.97x 10 8 


2.63xl0 12 


5.59xl0 14 


IIB 


;s,4; 


-8.10x 10 6 


2.14xl0 10 


5.57xl0 14 


IIIA 


'6,2^ 


3.93x 10 11 


-4.09x 10 11 


2.87x 10 14 


IIIA 


'8,4^ 


3.19x 10 9 


-3.22x 10 9 


2.85x 10 14 


IIIB 


'6,2^ 


3.85x 10 11 


-3.99x 10 11 


2.99x 10 14 


IIIB 


;s,4; 


3.13x 10 9 


-3.25x 10 9 


2.97x 10 14 



Table 2: Heavy right-handed Majorana neutrino masses Mj for degenerate models 
(IA,IB,IC), inverted models (IIA, IIB) and normal hierarchical models (IIIA, IIIB), with 
tan 2 012=0.5, using neutrino mass matrices given in Appendix A. The entry (m, n) in ttilr, 
indicates the type of Dirac neutrino mass matrix taken as charged lepton mass matrix (6,2) 
or up quark mass matrix (8,4), as explained in the text. 



rfhjGeV) 
1.19xl0" 9 



~Wh\! 



Type 



m,n 



4.78x IO" 6 



fei 



ei 



Y, 



B 



IA 
IA 



6,2 

8,4 



1.19xl(T 9 3.87x 10- 8 



9.3xl0- 5 
9.3xl0- 5 



1.53XKT 7 
4.14X10- 9 



1.55xl0" 13 
4.16xl0- 15 



5.31x IO" 7 
4.30x 1Q- 9 



4.46xl0" 16 
3.62xlQ- 18 



IB 
IB 



6,2 

8,4 



3.97x 
3.97x 



10 -io 

io- 10 



2.83xl0- 4 
2.83xl0- 4 



1.36xl0" 21 

l.lOxlO- 23 



5.31x IO" 7 
4.30x IO" 9 



2.49xl0" 15 
2.16xl0- 16 



7.62xl0- 21 
6.62xl0- 22 



IC 
IC 



4,2 
8,4 



3.97x 
3.97x 



10 -io 

io- 10 



2.83xl0- 4 
2.83xlQ- 4 



io- 11 



5.31x IO" 7 
4.30x IO" 9 



2.95xl0" 3 
2.95x10-3 



1.56xl0" 12 
1.26xl0- 14 



4.98xlO- 1Y 
4.04xl0- 19 



IIA 
IIA 



6,2 

8,4 



4.95x 

4.95x 



5.01x IO" 6 
4.06x IO" 8 



2.69xl0" 7 
2.18xl0- 9 



2.57xl0" 12 
2.07xl0- 14 



IIB 
IIB 



6,2 

8,4 



1.08x 
1.52x 



io- 12 
io- 10 



8.83xl0- 4 
8.83X10- 4 



7.51xl0- 3 
6.13xl0- 5 



4.59x IO" 3 
3.62x IO" 5 



9.27x IO" 9 
7.28x IO" 11 



IIIA 
IIIA 



6,2 

8,4 



5.80x 
5.80x 



io- 10 
io- 10 



1.82xl0- 4 
1.82xl0- 4 



7.51x10-3 
6.13xl0- 5 



4.91 x IO" 3 
3.88x IO" 5 



9.66x IO" 9 
7.59x IO" 11 



IIIB 
IIIB 



6,2 

8,4 



5.93x 
5.93x 



io- 10 
io- 10 



1.83xl0- 4 
1.83xl0- 4 



Table 3: Values of CP asymmetry and the baryon asymmetry for degenerate models 
(IA, IB, IC), inverted hierarchical models (IIA, IIB) and normal hierarchical models (IIIA, 
IIIB) with for tan 2 9\2 =0.50, using mass matrices given in Appendix A. The entry (m, n) 
indicates the type of Dirac mass matrix as explained in the text. 



Type (m,n) T% in <T R < T% ax (GeV) Mf in < M 7 < Mf ax (GeV) 


IA (6,2) 2.61 x 10 5 <T R < 1.22 x 10 6 2.44 x 10 8 < M 7 < 1.15 x 10 9 
IA (8,4) 7.80 x 10 4 <T R < 9.86 x 10 3 1.97 x 10 6 < M 7 < 2.49 x 10 5 


IB (6,2) 2.97 x 10 13 <T R < 4.05 x 10 5 9.10 x 10 7 < M 7 < 0.49 

IB (8,4) 2.97 x 10 13 < T R < 3.28 x 10 3 6.56 x 10 5 < M 7 < 7.20 x 10~ 5 


IC (6,2) 5.29 x 10 i3 <T R < 4.05 x 10 5 8.05 x 10 7 < M 7 < 2.74 

IC (8,4) 4.97 x 10 11 <T R < 3.28 x 10 3 6.56 x 10 5 < Mj < 4.35 x 10~ 3 


IIA (6,2) 6.80 x 10 1U <T R < 3.25 x 10 6 6.50 x 10 8 < M 7 < 3.09 x 10 4 
IIA (8,4) 6.80 x 10 10 <T R < 2.64 x 10 4 5.26 x 10 6 < M 7 < 2.03 


IIB (6,2) 1.22 x 10 6 <T R < 9.99 x 10 6 1.99 x 10 9 < M 7 < 1.64 x 10 10 
IIB (8,4) 1.22 x 10 6 <T R < 8.1 x 10 4 1.62 x 10 7 < M 7 < 1.08 x 10 6 


IIIA (6,2) 2.80 x 10 5 <T R < 3.93 x 10 9 7.86 x 10 11 < M 7 < 1.10 x 10 17 
IIIA (8,4) 2.88 x 10 5 <T R < 3.19 x 10 7 6.38 x 10 9 < M 7 < 7.06 x 10 12 


IIIB (6,2) 2.57 x 10 4 < T R < 3.85 x 10 9 7.70 x 10 11 < M 7 < 1.15 x 10 17 
IIIB (8,4) 2.64 x 10 4 < T R < 3.13 x 10 7 6.26 x 10 9 < M 7 < 7.42 x 10 12 



Table 4: Theoretical bounds on reheating temperature Tr and inflaton mass M 7 in 
non-thermal leptogenesis, for all neutrino mass models described in Tables 1-3. 



3 Non-thermal leptogenesis 

We next apply the neutrino mass models discussed in section 2 (Tables 1-3) 
to non-thermal leptogenesis scenario [7] where the right-handed neutrinos are 
produced through the direct non-thermal decay of the inflaton. We follow 
the standard procedure oulined in ref. [20] where non-thermal leptogenesis 
and baryon asymmetry in the universe had been studied in different neutrino 
mass models whereby some mass models were excluded using bounds from 
below and from above on the inflation mass and reheating temperature after 
inflation. Though we adopt similar analysis, the texture of the neutrino 
mass models considered here are different and hence the conclusions are also 
expected to be different. 

We start with the inflation decay rate given by 



IV = T(0 - NiNy 



I A; 



4tt 



-Mr 



(10) 



where \ are the Yukawa coupling constants for the interaction of three heavy 
right-handed neutrinos iV, with the inflaton of mass M 7 . The reheating 



temperature after inflation is given by the expression, 

t r =(-^-) 1/ \t,m p )^ (ii) 

\ 27r 9* J 

where Mp ~ 2.4 x 10 18 GeV is the reduced Planck mass[21] and g* is the 
effective number of relativistioc degrees of freedom at reheating temperature. 
For SM we have #* = 106.75 and for MSSM g* = 228.75. If the inflaton 
dominantly couples to N, the branching ratio of this decay process is taken 
as BR ~ 1, and the produced baryon asymmetry of the universe can be 
calculated by the following relation [22], 

Ys = ^ = CY L = C \^ (!2) 

where Yl is the lepton asymmetry generated by CP- violating out-of-equilibrium 
decays of heavy neutrino JVi and Tr is the reheating temperature. The frac- 
tion C has the value C = -28/79 for SM and C = -8/15 in the MSSM. 

The above expression (12) of the baryon asymmetry is supplemented 
by two more boundary conditions [20]: (i) lower bound on inflaton mass 
Mj > 2Mi coming from allowed kinematics of inflaton decay, and (ii) an 
upper bound for the reheating temperature Tr < O.OlMi coming from out- 
of-thermal equilibrium decay of JVi. Using the observed central value[2] of 
the baryon asymmetry Yb — — — 8.7 x 10~ n and theoretical prediction of 
CP asymmetry e in Table 3, in equation (12), one can establish the relation 
between Tr and Mi for each neutrino mass model. 

The right-handed neutrino mass M\ from Table 2 and the CP asymmetry 
e from Table 3 for all neutrino mass models, are used to calculate the bounds: 
T min < t r < T^ ax and M™ in < M 7 < M 7 max in Table 4 following eq.(12) 
along with other two boundary conditions cited above. Only those models 
which satisfy the constraint T™* > T R mn could survive in the non-thermal 
leptogenesis. These models are identified as IA with (6,2), IIB with (6,2), 
III (A,B) with (6,2) and III(A,B) with (8,4) where (m,n) refers to the type 
of Dirac neutrino mass matrix. From Table 4 it is seen that inflationary 
models in which Mj ~ 10 13 GeV, like e.g., chaotic or natural inflation, are 
compatible only with normal hierarchical model III (A, B) with (6,2). In fact 
with T R = 10 6 GeV, we get M, = 2.8 x 10 13 GeV., T^ = 2.85 x 10~ 6 GeV, 
and |Ai| = 1.13 x 10 -8 which are compatible with chaotic inflationary model. 

In supersymmetric models, the gravitino problem [8] can be avoided pro- 
vided that the reheating temperature after inflation is bounded from above 

10 



in a certain way, namely Tr < (10 6 — 10 7 ) GeV. In fact the reheating temper- 
ature Tr = 10 6 GeV is relevant in order to realise the weak scale gravitino 
mass rri3/2 ~ 100 GeV without causing the gravitino problem. Even this 
reheating temperature is relaxed for two order Tr = 10 7 GeV, we would 
have Mj ~ 10 11 GeV in normal hierarchy type III(A,B) with (8,4). We con- 
clude that the only surviving model in this analysis is the normal hierarchical 
model (III). 

4 Summary and discussion 

To summarise, we first parametrise the light left-handed Majorana neutrino 
mass matrices describing the possible patterns of neutrino masses viz, degen- 
erate, inverted hierarchical and normal hierarchical, which obey the \x — r 
symmetry having tribimaximal mixings. As a first test these mass matrices 
predict the neutrino mass parameters and mixings consistent with data, and 
all the input parameters are fixed at this stage. In the next stage these mass 
matrices are employed to estimate the baryon asymmetry in both thermal as 
well as non-thermal leptogenesis scenario. We use the CP violating Majorana 
phases derived from right-handed Majorana mass matrix and two possible 
forms of Dirac neutrino mass matrices as either charged lepton mass ma- 
trix or up-quark mass matrix in the calculation. The overall analysis shows 
that normal hierarchical model appears to be the most favourable choice 
in nature. The present analysis though phenomenological may serve as an 
additional criteria to discard some of the presently available neutrino mass 
models and neutrino mass ordering patterns. There are some suggestions in 
the literature [24] for inverted hierarchical model to enhance the estimation 
of baryon asymmetry if m 3 is increased. The present investigation has taken 
care of the maximum allowed non-zero value of m^ ~ 0.033 eV in case of in- 
verted hierarchy type IIB model. Our result also differs from a recent study 
in nonthermal leptogenesis with strongly hierarchical right-handed neutri- 
nos [25] where the mass of the lightest right handed neutrino Mi < 10 6 GeV. 
There are some propositions[26] for probing the reheating temperature at 
the Large Hadron Collider and this hopefully decides the validity of thermal 
leptogenesis. 



11 



Appendix A 

Classification: 

We first list here for ready reference to the classification of neutrino mass 
models, the zeroth-order left-handed Majorana neutrino mass matrices with 
texture zeros, itill, corresponding to three models of neutrinos given in the 
text, viz., degenerate (Type [I]), inverted hierarchical (Type [II]) and normal 
hierarchical (Type [III]). 



Type 



m LL 



m 



diag 
LL 



[IA] 



I ° j- 2 £\ 

1 1 _1 

V2 2 2 

, J- -I I / 
\ V2 2 2/ 



itlq Diag(l,—l,l)mo 



[IB] 



[IC] 




Diag(l, 1, l)mo 



Diag(l, 1, — 1)toq 



[IIA] 



[IIB] 




Diag(l, 1, 0)m 



Diag(l, — 1, 0)mo 



[III] 




Diag(0, 0, l)mo 
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Parametrisation with two parameters for tribimaximal 
mixings: 

Left-handed Majorana neutrino mass matrices which obey fi — r symme- 
try[10,23] have the following form 



m LL 



This predicts an arbitrary solar mixing angle tan 29\2 = | (x-z-w) I ' wm ^ e the 
predictions on atmospheric mixing angle is maximal (#23 = 7r /4) and Chooz 
angle zero. We parametrise the mass matrices (with only two parameters) 
whereby the solar mixing is fixed at tribimaximal mixings for all possible 
patterns of neutrino mass models: 

l.Deg Type A [lA](mj = mi, —m2,m 3 ) 

/5 1 -25 2 -5 1 -5 1 \ 

m LL =\ -81 \~8 2 -\-8 2 \m 
V -*i -I-S2 \~8 2 ) 

with input values: <5i=0.66115, <5 2 =0.16535,m o = OAeV. 

2.Deg Type B [IB](toj = mi,m 2 ,m 3 ) 

1 — 81 — 28 2 8\ 8\ 

m LL = jj 81 1- 8 2 -8 2 I m c 

81 —82 1 — 82 

with input values: <5i=8.314xl(r 5 ,<52=0.00395,m o =0.4eV. 

3.Deg Type C [IC](mj = m 1 ,m2, —m 3 ) 

1 — 81 — 282 8\ 8\ 

m LL = J 81 -8 2 l-8 2 \m c 

81 1 — 82 —82 

with input values: 5i=8.314xl0" 5 ,<5 2 =0.00395,m o =0.4eV. 
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4: Inverted Hierarchical mass matrix with m 3 7^ 0: 

1 - 2e -e -e 

m LL (IH)= I -e 1/2 1/2 - 77 | m . 

-e I/2-77 V 2 



Inverted Hierarchy with even CP parity in the first two mass eigenvalues 
[IIA] (mi = mi,m 2 ,m 3 ): 7//e=1.0,7/=0.0048,m = 0.05eV. 
Inverted Hierarchy with odd CP parity in the first two mass eigenvalues [IIB] 
(rrii = mi,-m 2 ,m 3 ): ?y/e=1.0,77=0.6607,mo = 0.05e^. 

5: . Normal Hierarchical mass matrix Case (i) with m(l, 1) 7^ type- 
[IIIA]: 




m LL {NH) 



with input values: 7//e=0.0,e=0.175,m = 0.029CV. 

6: Normal Hierarchical mass matrix Case (ii) with m(l, 1) = type- 
[IIIB]: 

/ -c -c \ 

m LL (NH)=\ -e 1-e — 1 + 77 m 
\ — e —I+77 1-e / 

with input values: rj/e—0. 0, e=0.164,m = 0.028CV. 

The textures of mass matrices for inverted hierarchy (IIA, IIB) as well as 
normal hierarchy (IIIA, IIIB) have the potential to decrease the solar mixing 
angle from the tribimaximal value, without sacrificing fi — r symmetry. This 
is possible through the identification of 'flavour twister' 77/e 7^ [23]. 
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